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Courbes  et  Surfaces  - Ellipse  et  Ellipsoïde 
de  révolution 

1.1  Courbes  et  Surfaces 

1.1.1  Courbes  Gauches 

* Le  trièdre  de  Frenêt. 

* Définition  de  la  courbure. 

* Calcul  dans  le  cas  d’une  courbe  plane. 

1.1.2  Surfaces  paramétrées 

* Plan  tangent,  vecteur  unitaire  normal. 

* 1ère  forme  quadratique  fondamentale,  courbes  orthogonales  symétriques. 

* Courbe  tracée  sur  une  surface,  trièdre  de  Darboux. 

* Théorème  de  Meusnier. 

* Directions  et  courbures  principales,  courbure  d’une  section  normale. 

* Courbures  principales  d’une  surface  de  révolution. 

1.2  L’Ellipse  et  l’ellipsoïde  de  révolution 

* Définition. 

* Définitions  concernant  l’aplatissement,  les  excentricités. 

* Définitions  des  latitudes  : 

- paramétrique  : if), 
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- géographique  : <p, 

- géocentrique  : u, 

1.2.1  Exercices 

1.  Calculer  les  composantes  du  vecteur  normal  extérieur  à l’ellipsoide,  en  déduire  les 
relations  (dans  les  deux  sens)  entre  les  lignes  trigonométriques  de  ip  et  celles  de  i/j. 

2.  Donner  les  équations  paramétriques  de  l’ellipse  et  de  l’ellipsoide  en  fonction,  re- 
pectivement,  de  cp  et  de  À et  <p. 

3.  Etablir  une  relation  différentielle  entre  ÿ et  p>. 

4.  Calculer  la  différentielle  d/3  de  l’arc  d’ellipse  en  fonction  de  <p,  puis  la  première 
forme  quadratique  de  l’ellipsoide. 

5.  Calculer  les  courbures  principales  de  l’ellipsoide  de  révolution. 

6.  Trouver  la  coordonnée  curviligne  de  l’ellipsoide  de  révolution  qui  forme  avec  la 
longitude  un  couple  de  coordonnées  symétriques  et  qui  s’annulle  le  long  de  l’équateur. 


1.3  Calcul  d’un  arc  d’ellipse 


1.3.1  Intégrales  de  Wallis 

Soit  à calculer  : 


On  pose  : 


pÇl 

W2 p = / sin2pijjduj 

J 0 

pÇl 

Ip- 2(H)  = / sinp~2cjjœs2ujduj 

J 0 


1.  Etablir  les  formules  suivantes  : 


Wp  = 


(M 

1 

1 

CM 

1 

£ 

II 

£ 

(1.3.1) 

l)/p_2  = sinp  1Dco-sD  + Wp 

(1.3.2) 

P~1Wp-2  1 sinp~1ü,cosQ 

(1.3.3) 

P P 

2.  Préciser  la  valeur  de  TT q • et  proposer  un  programme  (en  Matlab)  de  calcul  de  W 2P ■ 


1.3.2  Calcul  d’un  arc  d’ellipse  méridienne 


1.3.  CALCUL  D’UN  ARC  D’ELLIPSE 
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avec  : 


o(l  6 ) 2 -,  2-2 

p = tj , w = 1 — e sm  <p 


w d 


1.  Développer  te-3  suivant  les  puissances  croissantes  de  esirnp. 

2.  Calculer  {3((p)  en  fonction  des  W2P(<p). 

3.  Majorer  l’erreur  de  calcul,  lorsqu’on  arrête  le  développement  au  terme  e2n.  Calculer 
n si  l’on  recherche  la  précision  du  millimètre  sur  f3,  quelle  que  soit  la  latitude  cp  entre 


-f  et  f • 

4.  Proposer  un  organigramme  de  calcul. 

5.  Envisager  la  solution  du  problème  inverse  : calcul  de  ip  connaissant  (3. 
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Chapitre 

Passage  des  coordonnées  cartésiennes  en 
ellipsoidiques 

2.1  Introduction 

On  considère  un  référentiel  géodésique  ( 0,0X,0Y,0Z ) avec  un  ellipsoide  de  ré- 
férence £ de  paramètres  a,  e2  où  respectivement  a le  demi-grand  axe  et  le  carré  de  la 
première  excentricité.  Les  coordonnées  cartésiennes  ( X , Y,  Z)  d’un  point  M((p,  A,  h)  sont 
données  par  : 

X = ( N + h)cosipcos\ 

M < Y = (N  + h)cos<psin\  (2.1.1) 

Z = (N(  1 — e2)  + h)sinip 

On  veut  étudier  le  passage  de  (X,  Y,  Z)  à (cp,  A,  h). 

1.  Montrer  en  considérant  que  Z > 0,  h vérifie  : 

h > —(N(l  — e2)  (2.1.2) 

Le  calcul  de  A est  facile,  et  on  a : 

Le  calcul  de  tp,  h est  plus  complexe 
suivantes  : 

a - les  algorithmes  itératifs, 
b - les  algorithmes  finis, 
c - les  développements  limités. 
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t9(  A)  = ^ (2.1.3) 

. Sa  résolution  peut  se  faire  par  les  trois  méthodes 


UCHAPITRE  2.  PASSAGE  DES  COORDONNÉES  CARTÉSIENNES  EN  ELLIPSOIDIQUES 

2.1.1  Les  Algorithmes  Itératifs 

Les  algorithmes  itératifs  conduisent  à résoudre  une  équation  de  la  forme  : 


x = f(x) 


(2.1.4) 


où  x est  une  inconnue  auxiliaire.  Soit  x une  solution  de  (2.1.4).  Partant  d’une  solution 
approchée  xq  de  x,  on  calcule  successivement  : 


xi  = f(x  0) 
X2  = f(x  i) 


Xi  = f(Xi-l) 


La  méthode  converge  si,  pour  un  ensemble  de  voisinages  V)  de  la  solution  x : 


Vi  = f(Vi- 1)  C Vi-! 


(2.1.5) 


On  dit  que  / est  une  fonction  contractante  au  voisinage  de  x. 

Si  / est  continue  et  dérivable  au  voisinage  de  x,  la  condition  de  convergence  est  : 


3 un  voisinage  V de  x 
3k  G [0, 1] 


tels  que  x G V =>-  |//(æ)|  < k 


(2.1.6) 


En  d’autres  termes,  la  méthode  itérative  converge  si  \ f' (x)\  est  majorée  par  un  coefficient 
k inférieur  à 1 au  voisinage  de  la  solution.  La  divergence  de  la  méthode  n’implique  pas 
l’inexistance  de  la  solution. 


Etude  de  la  précision  de  la  Méthode 

de  : 


xi  = f(x  0) 
X2  = f(x  l) 


Xi  = f(Xi- 1) 


2.1.  INTRODUCTION 

on  tire  : 

x-xi=  f(x)  - f{x0)  = (x-  xo)f(0i) 
x-x2  = f(x)  - f(x i)  = {x-  Xl)f'(02) 

X - Xi- 1 = f(x)  - f(Xi-2 ) = (x  - Xi_2)//(6»j_l) 

X-Xi  = f(x)  - f(Xi- 1)  = (x  - Xi-l)f'(0i) 

avec,  à chaque  fois  9i  G [xj,x].  On  en  déduit  les  i + 1 inégalités  : 

\x  — x0|  < \x  — x0| 

|x-xi|  < k\x  — xo| 

\X  — Xi- 1|  < k\x  — Xi— 2 1 
\x-Xi\  < k\x  — Xi- 1| 

Soit  par  multiplication  membre  à membre  : 

\x-Xi\  < kl\x-  x0\  (2.1.8) 

Comme  la  fonction  est  contractante,  on  peut,  dans  cette  équation,  remplacer  \x  — xq\  par 
\x'0  — Xo|,  x'Q  étant  une  valeur  approchée  de  x formant,  avec  xq  un  encadrement  de  x : 

X G [x0,x'o] 

L’erreur  de  la  méthode  est  alors  majorée  par  : 

\x  — Xi\  < kl \x'0  — xo|  (2.1.9) 

Le  nombre  i d’itérations  nécessaires  pour  calculer  x avec  une  précision  ex  donnée  à 
l’avance  est  obtenu  par  : kl\x'0  — xq\  < ex  soit,  en  se  souvenant  que  k < 1 

(2.1.10) 

Etapes  de  l’Analyse  et  la  programmation  d’une  Méthode  itérative  : 

L’Analyse  et  la  programmation  d’une  Méthode  itérative  peuvent  comporter  les  élé- 
ments suivants  : 

Démonstration  de  l’équation  employée  : x = f(x). 
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(2.1.7) 
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- Calcul  et  majoration  de  /',  eventuellement  étude  des  cas  où  \ f'(x)\  > 1. 

- Calcul  de  i,  nombre  d’itérations  nécessaires.  Ce  nombre  peut  être  fixé  une  fois  pour 
toutes,  ou  calculé  par  le  programme. 

- Calcul  de  ip  et  de  h par  une  formule  peu  sensible  aux  erreurs  sur  x. 

- Rédaction  d’un  organigramme  et  d’un  programme  en  Matlab,  accompagnés  d’une 
notice  d’emploi  sur  les  limites  d’emploi,  le  temps  d’éxecution  et  l’encombrement 
en  machines. 

- Un  jeu  d’essai  au  moins. 


Les  Méthodes  itératives  proposées  : 

1.  tp  = Arctg  (-  + Ne2sin(f\  avec  p = ^X2  + Y2. 


2.  ip  = Arctg 


\P 

Z 

P 


P ) 


1 - 


Ne2cos(p\ 


-r 


P 


avec  p = VX2  + Y2. 
Z 


Ne2sin2ip  a — r 

3.  (p  = ip  + Arcsm , avec  w = Arctg—,  =,  r = \/ pz  + Z . 

2 r yjx2  + Y2 

2.1.2  Les  Algorithmes  finis 

Les  algorithmes  finis  conduisent  à la  résolution  d’une  équation  du  4ème  degré  : 

X4  + CllX3  + CL2X2  + CI3X  + CL4  = 0 (2.1.11) 

où  x est  une  variable  auxiliaire.  O11  donne  ci-dessous  quelques  indications  pour  la  réso- 


lution de  l’équation  (2.1.11). 

On  élimine  le  terme  du  troisième  degré  par  la  transformation  linéaire  : 

. «i 
V = x + — 


(2.1.12) 


L’équation  (2.1.11)  devient  : 


V4  + a'2y2  + a'3y  + 0,4  = 0 


(2.1.13) 


On  abaisse  ensuite  le  degré  de  cette  équation  en  posant  : 

2 y = u + v + w (2.1.14) 

Entre  les  variables  indépendantes  u,v,w,  on  peut  encore  imposer  deux  relations,  par 
exemple  : 

u 2 + v2  + w2  = — 2a'2 
u.v.w  = —a3 


(2.1.15) 

(2.1.16) 


2.1.  INTRODUCTION 
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En  utilisant  les  fonctions  symétriques  des  racines,  on  montre  que  a2,  v2,w2  sont  solutions 
de  l’équation  : 

z3  + 2 a2z2  + (a22  - 4 a'A)z  - a'3  = 0 (2.1.17) 


On  fait  disparaître  le  terme  du  second  degré  par  le  changement  : 

. 2a'2 

et  on  résout  l’équation  par  la  méthode  de  Cardan  (Chercher  dans  le  Web). 


(2.1.18) 


L’équation  (2.1.17)  admet  trois  solutions  dans  le  corps  des  nombres  complexes  C : 
u2,v2,w2.  Ce-ci  fournit  à priori  8 possibilités  de  calculer  y par  (2.1.14).  Cependant,  le 


signe  du  produit  uvw  est  imposé  par  (2.1.16),  et  l’équation  (2.1.11)  n’admet  bien  que 
quatre  solutions  parmis  lesquelles  il  faudra  choisir. 


L’Analyse  et  la  Programmation  d’une  Méthode  Finie  : 

L’Analyse  et  la  programmation  d’une  méthode  finie  comportera  les  éléments  sui- 
vants : 

Démonstration  de  l’équation  employée. 

Le  calcul  des  racines  et  la  discussion. 

- Le  calcul  de  tp  et  de  h. 

La  rédaction  d’un  organigramme  et  d’un  programme  en  Matlab,  accompagnés 
d’une  notice  d’emploi  sur  les  limites  d’emploi,  le  temps  d’éxecution  et  l’encom- 
brement en  machines. 

- Un  jeu  d’essai  au  moins. 


Les  Méthodes  Finies  Proposées  : 


1. 


2. 


On  écrit  que,  simultanément,  le  point  m appartient  à l’ellipse  méridienne  et  que  sa 

distance  à M est  extrémale.  C’est  un  problème  d’extrémum  lié.  On  se  ramènera  à 

l’équation  (2.1.11)  dans  laquelle  x est  le  multiplicateur  de  Lagrange  (Fig|2.1.T). 

a 2 

On  écrit  que  la  pente  de  mM  et  celle  de  ml  valent  tgp>  = tgu  où  u est  la  latitude 

0“ 

géocentrique  du  point  m.  Il  vient  : 


Z — Zq  a2  Zq 
R — Rq  b2  Rq 


(2.1.19) 


On  élimine  Zq  à l’aide  de  l’équation  de  l’ellipse  et  on  obtient  une  équation  du  4ème 
degré  en  Rq  (Fig.  |2.1.2  ). 
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Figure  2.1.1  - La  Méthode  Finie  1. 


Z 


Figure  2.1.2  - La  Méthode  Finie  2. 


2.1.  INTRODUCTION 
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2.1.3  Méthode  des  Développements  limités 

x étant  une  inconnue  auxiliaire,  et  t un  paramètre  facilement  calculable  en  fonction 
des  données,  on  cherche  les  ai  tels  que  : 

i=n 

x = ^jaitn  (2.1.20) 

i= 0 

Les  ai  sont  des  coefficiens  fonctions  des  données  X , Y , Z , a e te2;  ils  sont  connus  jusqu’à 
un  certain  rang,  ou  calculables  par  un  programme. 

Il  peut  être  difficile  de  majorer  l’erreur  de  la  méthode,  suivant  les  cas.  Le  développe- 
ment peut  aussi  diverger  pour  certaines  valeurs  de  t. 

L’Analyse  et  la  Programmation  : 

L’Analyse  et  la  programmation  d’une  méthode  par  développement  limité  peut  com- 
porter : 

1.  La  démonstration  des  formules  employées. 

2.  Une  étude  analytique  de  précision  ou,  si  impossible,  une  comparaison  avec  d’autres 
méthodes. 

3.  La  rédaction  d’un  organigramme  et  d’un  programme  en  Matlab,  accompagnés 
d’une  notice  d’emploi  sur  les  limites  d’emploi,  le  temps  d’éxecution  et  l’encom- 
brement en  machines. 

4.  Un  jeu  d’essai  au  moins. 

Le  Développement  proposé  : 

On  pose  : 


R = (N  + h)cos<p 

Z = (IV(1  — e2)  + h)simp 
R 

x = zUj'-p 

v2  = {1-  e2)Z2 
+ _ e2 

~~  VR2  + v'2 

p = \/  R2  + x2u2 
R2 

R2  + u2 


(2.1.21) 

(2.1.22) 

(2.1.23) 

(2.1.24) 

(2.1.25) 

(2.1.26) 
(2.1.27) 


c = 
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1.  Montrer  que  : 


2.  Montrer  que  x vérifie  l’équation 


9 ax 

x — 1 = e2 — 


(x  - 1 )2(R2  + v2x2)  = eAa2x2 


(2.1.28) 


(2.1.29) 


On  pose  alors  : 

i=n 

x = x = ^2aitn  (2.1.30) 

i= 0 

3.  Trouver  les  ai  pour  0 < i < 2 par  identification  des  deux  membres  de  l’équation  en  x. 
On  admettra  ensuite  : 


«3  = 


04  = 2c  — 9c2  + 8c3 


(2.1.31) 

(2.1.32) 


Deuxième  partie 


Exercices  et  Problèmes 
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Chapitre 

Exercices  et  Problèmes 


3.1  Trigonométrie  Sphérique 

Exercice  n°  1 : Calculer  l’azimut  d’une  étoile  de  déclinaison  ô = +5°  quand  sa  dis- 
tance zénithale  est  de  80°  pour  un  observateur  situé  à la  latitude  cp  = 56°. 

Exercice  n° 2 : En  appliquant  au  triangle  de  position  les  formules  de  trigonométrie 
sphérique  montrer  que  l’on  peut  calculer  l’angle  horaire  AHC  du  coucher  d’un  astre  par  : 
cosAHc  = —tgtp.tgô. 

Exercice  n° 3 : Soit  un  triangle  sphérique  ABC.  On  donne  les  éléments  suivants  : 

- Â = 80.1643  3 gr, 

- è = 55.7735  1 gr, 

- C = 64.0626  1 gr, 

- AC  = 20.1357  km, 

- AB  = 22.1435  km. 

1.  Calculer  a = Â + È + C. 

2.  Déterminer  e l’excès  sphérique  de  ce  triangle. 

3.  Calculer  la  fermeture  du  triangle  ABC,  donnée  par  : 

f = a — 200.0000  Ogr  — e 


Exercice  n°4  : Soit  (S2)  une  sphère  de  rayon  égal  à 1.  Soit  un  carré  sphérique 
ABCD  de  côté  a (arc  de  grand  cercle).  On  note  a = A = B = C = D. 
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CHAPITRE  3.  EXERCICES  ET  PROBLEMES 


Figure  3.1.1  - Les  coordonnées  de  Cassini-Soldner 


1.  Montrer  que  : 


cosa  = cotg1  2 


a 

2 


2.  Donner  l’expression  de  la  diagonale  d = l'arc  AC. 


Problème  n°l  : Soit  (S2)  une  sphère  de  rayon  égal  à 1 et  de  centre  le  point  D. 
Un  point  M de  (S2)  a pour  coordonnées  (<p,-\).  On  appelle  les  coordonnées  de  Cassini- 
Soldner T 2 de  M les  angles  (Fig|3.1.ï|  : 


- L = no,  OH, 

- h = nifjüi. 

1.  Déterminer  les  relations  liant  L,H  à (p,  À. 

2.  Inversement,  donner  les  relations  liant  ip,  A à L,H. 

Problème  n° 2 : Au  lieu  M de  latitude  ip  = 38°  Nord,  on  observe  l’étoile  polaire  A 
de  déclinaison  ô = +89°  et  d’ascension  droite  a — +2  h 13  mn  52.90  s. 

1.  Donner  sur  un  graphique,  les  éléments  du  triangle  sphérique  P AAI  où  P est  le 
pôle  Nord. 

2.  Sachant  que  l’heure  sidérale  locale  H SL  est  égale  au  moment  de  l’observation  à 
6 h 37  mn  19.72  s,  calculer  l’angle  horaire  AH. 

3.  En  appliquant  la  formule  des  cotangentes,  montrer  que  l’azimut  Az  de  l’étoile  est 
donné  par  la  formule  : 


tgAz 


sinAH 

cosAHsinip  — cos<ptgô 


1.  César-François  Cassini  (1714-1784)  : Astronome  et  géodésien  français. 

2.  Dr  Johann  Georg  von  Soldner  (1776-1833)  : Mathématicien  et  astronome  bavarois. 
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4.  Calculer  alors  l’azimut  Az. 

5.  Calculer  la  distance  zénithale  z de  l’étoile. 
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3.2  Astronomie  de  Position 

Exercice  n°  1 : Au  lieu  de  latitude  <p  = 36°54/  Nord,  on  veut  calculer  les  hauteurs  h \ 
et  li-2  de  l’étoile  polaire  de  déclinaison  ô = +89°  respectivement  à son  passage  supérieur 
et  à son  passage  inférieur  au  méridien  du  lieu.  Déterminer  h\  et  h^. 

Problème  n°  1 : 1.  En  un  lieu  de  latitude  ip  quelles  sont  les  étoiles  : 

- qui  ne  se  couchent  pas  ( qui  sont  toujours  visibles), 

- qui  ne  sont  jamais  visibles. 

Traiter  le  cas  : lieu  dans  l’hémisphère  nord. 

2.  Quelle  est  la  condition  pour  qu’une  étoile  culmine  au  zénith? 

3.  Cas  particulier  du  soleil  : la  déclinaison  du  soleil  varie  de  — 23°27/  à +23°27/  au 
cours  de  l’année.  On  appelle  jour  le  moment  pendant  lequel  le  soleil  est  au-dessus  de  l’ho- 
rizon, nuit  lorsque  le  soleil  est  au-dessous  de  l’horizon,  midi  l’instant  de  la  culmination, 
minuit  l’instant  du  passage  inférieur. 

a)  Montrer  qu’au  moment  des  équinoxes  le  jour  et  la  nuit  sont  d’égale  durée  quel  que 
soit  le  lieu. 

b)  Montrer  qu’à  l’équateur,  quelle  que  soit  la  date  le  jour  et  la  nuit  sont  d’égale  durée. 

c)  Au  moment  du  solstice  d’hiver  quels  sont  les  lieux  : 

- où  il  fait  constamment  jour, 

- où  il  fait  constamment  nuit. 

Mêmes  questions  au  moment  du  solstice  d’été. 

d)  Quels  sont  les  lieux  de  la  terre  où  le  soleil  culmine  au  zénith  au  moment  du  solstice 
d’hiver.  Même  question  au  moment  du  solstice  d’été. 

e)  Quels  sont  les  lieux  de  la  terre  où  au  moins  une  fois  dans  l’année  le  soleil  culmine 
au  zénith. 

Problème  n° 2 : Une  station  astronomique  est  située  en  un  lieu  de  coordonnées 
géographiques  : <p  = +45°  00';  A = +7  h 20  mn. 

En  ce  lieu,  on  observe  une  étoile  A de  coordonnées  équatoriales  : 

a = +11  h 13  mn;  ô = 30°  00'. 

L’observation  se  fait  le  jour  de  l’équinoxe  de  printemps  le  21  mars  à 0 heure  TU. 
L’heure  sidérale  de  Greenwich  est  11  h52mn. 

1.  Calculer  l’heure  sidérale  locale  du  lever  et  du  coucher  de  l’étoile  A au  lieu  considéré. 

2.  En  déduire  l’heure  TU  du  lever  et  du  coucher  de  l’étoile  au  heu  considéré. 
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Remarque  : on  choisira  le  coucher  qui  a lieu  après  le  lever. 

Problème  n° 3 : En  un  lieu  de  latitude  43°,  521  et  de  longitude  +0  h 20  mn  57  s,  on 
cherche  à pointer  la  galaxie  d’Andromède  de  coordonnées  équatoriales  a = 0 h 40  mn,  5 = 
41°  00'  le  31  juillet  1992  à 21  hTU. 

On  donne  l’heure  sidérale  de  Greenwich  à 0 hTU  le  31/07/1992  : HSGohTU  = 
20  h 35  mn  28  s. 

1.  Calculer  l’heure  sidérale  locale  à 21  hTU. 

2.  En  déduire  l’angle  horaire  de  la  galaxie. 

3.  Calculer  la  distance  zénithale  de  la  galaxie  à 21  hTU. 

4.  Calculer  son  azimut  à cette  même  heure. 

Problème  n°4  : En  un  lieu  de  l’hémisphère  Nord  de  latitude  ip,  on  mesure  la  longueur 
de  l’ombre  portée  HC,  à midi  vrai  (passage  du  soleil  au  méridien),  par  une  tige  verticale 
HA  dont  l’extrémité  H est  sur  le  sol  supposé  horizontal. 

1.  Donner  l’expression  HC  en  fonction  de  HA  et  de  la  distance  zénithale  Dz  du 
soleil. 

2.  Donner  l’expression  de  HC  en  fonction  de  HA  et  de  p : 

- aux  équinoxes, 

- aux  solstices. 

3.  Quelle  doit  être  la  déclinaison  du  soleil  et  en  quels  lieux,  pourque  l’on  ait  HC  = 
HAÏ 

4.  En  un  lieu  de  latitude  c p = 47°  en  quelles  saisons  peut  on  avoir  HC  = HA. 

5.  Si  on  déplace  HB  le  long  d’un  méridien,  en  restant  dans  l’hémisphère  Nord,  existe- 
t-  il  au  cours  de  l’année  des  lieux  où  HC  = 0,  ou  HC  devient  infiniment  grand. 
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3.3  Courbes  et  Surfaces 


Exercice  n°  1 : Soit  l’hélice  circulaire  T paramétrée  par  : 

x = acost 
< y = asint 
z = bt 

où  a , b deux  constantes  positives. 

1.  Exprimer  les  composantes  des  vecteurs  T,  N,  B du  repère  de  Frénet. 

2.  Montrer  que  la  courbure  vaut  —5 — 3- . 

az  + bz 

3.  Montrer  que  la  torsion  vaut  — ^ 73 . 

a1  + b2 


Problème  n°  1 : Soit  la  courbe  ( C ) définie  par  les  formules  : 


M < 


x = at2 


y = at 3 

9 H 
Z=16at 


avec  a > 0 


1.  Calculer  l’abscisse  curviligne  s d’un  point  AI  quelconque  de  cette  courbe  lorsqu’on 
prend  pour  origine  des  arcs  l’origine  des  coordonnées  et  qu’on  prend  pour  sens  des  arcs 
croissants  celui  des  y croissants. 

2.  Déterminer  au  point  M les  vecteurs  unitaires  du  trièdre  de  Frenêt. 

3.  Calculer  le  rayon  de  courbure  et  les  coordonnées  du  centre  de  courbure. 

4.  Evaluer  la  torsion  en  M . 


Problème  n° 2 : Soit  (r)  la  surface  paramétrée  par  (u,v)  dans  M2  telle  que  : 


M (u,  v ) < 


X = u(  1 — u2)cosv 
Y = it(l  — u2)sinv 
Z = l-u2 


1.  Calculer  l’expression  de  ds2 . 

2.  Montrer  que  l’équation  cartésienne  de  (r)  est  : 


x2  + y2  = (1  - z)z 2 
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Exercice  n° 2 : Soit  la  surface  d’Enneper  : 


M(u,v ) < 


^ u 9 

X = u h uv 

§ 

V O 

y = v — — + vu 
Z = u1 2 3 4 5  — v2 


1.  Montrer  que  : 


ds2  = (1  + u2  + v2)2.(du2  + dv2) 


2.  Calculer  un  vecteur  unitaire  normal  à la  surface. 

3.  Montrer  que  la  surface  d’Enneper  est  de  courbure  moyenne  nulle  en  chaque  point. 


Exercice  n°3  : On  suppose  que  la  métrique  d’une  surface  donnée  est  : 

ds2  = A2  du2  + B2  dv2 , A = A(u,v),  B = B(u,v) 

1.  Montrer  alors  que  l’expression  de  la  courbure  totale  est  : 


' désigne  la  dérivation  partielle. 


Problème  n°3  : Soit  l’ellipse  (E)  définie  par  les  équations  paramétriques  : 


M < 


x = acosu 
y = bsinu 
avec  a > b > 0 


On  pose  : 

2 a2  — b2  ,2  a2  — b2 

1.  Calculer  la  position  sur  l’axe  des  abscisses  des  deux  points  F et  F'  appelés  foyers  tels 
que  MF  + MF'  = 2 a. 

2.  Montrer  que  le  produit  des  distances  des  foyers  à la  tangente  à l’ellipse  en  M est 
indépendant  de  u. 

3.  Donner  l’expression  de  ds. 

4.  Déterminer  les  expressions  des  vecteurs  unitaires  T et  Al  et  en  déduire  le  rayon 
de  coubure  de  l’ellipse. 

5.  Montrer  qu’il  passe  par  M deux  cercles  tangents  en  ce  point  à la  courbe  et  centrés 
sur  Ox,Oy  respectivement  (appelés  cercles  surosculateurs) . 
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6.  Que  deviennent  ces  cercles  lorsque  M est  un  sommet  de  l’ellipse. 


Problème  n°4  : On  définit  une  surface  (S)  par  les  équations  : 

X = u2  + v 


M(u,v ) < 


Y = u + v2 
Z = uv 


1.  Calculer  les  composantes  des  vecteurs  OMtu  et  OM'v. 

2.  Calculer  les  coefficients  E,F,G  de  la  première  forme  fondamentale  de  la  surface 
(5). 

3.  En  déduire  l’expression  de  ds2. 

4.  Les  coordonnées  (u,v)  sont-elles  orthogonales?  symétriques? 

5.  Calculer  un  vecteur  normal  de  (5). 


Problème  n° 5 : On  définit  une  surface  (E)  par  les  équations  : 


M(u,v ) < 


X = a.cosu.cosv 
Y = a.cosu.sinv 
Z = b.sinu 


avec  a,  b deux  constantes  positives. 

1.  Calculer  les  composantes  des  vecteurs  OM,u  et  OM,v. 

2.  Calculer  les  coefficients  E,F,G  de  la  première  forme  fondamentale  de  la  surface 

(E). 

3.  En  déduire  l’expression  de  ds2. 

4.  Les  coordonnées  (u,v)  sont-elles  orthogonales?  symétriques? 

5.  Calculer  un  vecteur  unitaire  normal  n de  (E). 

6.  Calculer  les  vecteurs  : 


OM"  OM"  OM" 

Ky±y±uu’>  Kyiyjuv’>  K-'iyjvv 


On  pose  : 


l = n.  om;iu, 


M = n.  OM'uv , N = n . OM'vv 


7.  Calculer  les  coefficients  L,  AI  et  N. 
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Problème  n° 6 : On  considère  la  surface  (T)  définie  par  les  équations  : 


M(u,  v ) < 


X = sinu.cosv 
Y = sinu.sinv 

U 

Z = cosu  + Logtg—  + ip(v) 


avec  est  une  fonction  définie  de  classe  C 1 de  v. 

1.  Donner  le  domaine  de  définition  de  la  surface  (T). 

2.  Montrer  que  les  courbes  coordonnées  v = constante  constituent  une  famille  de 
courbes  planes  de  (T)  et  que  leur  plan  coupe  (T)  sous  un  angle  constant. 

3.  Calculer  les  composantes  des  vecteurs  OMu  et  OMv. 

4.  Calculer  les  coefficients  E,  F,  G de  la  première  forme  fondamentale  de  la  surface 

(r). 

5.  En  déduire  l’expression  de  ds2. 

6.  Les  coordonnées  (u,v)  sont-elles  orthogonales?  symétriques? 

7.  On  suppose  pour  la  suite  que  îp(v)  = 0,  calculer  un  vecteur  unitaire  normal  n de 

r. 

8.  Calculer  les  vecteurs  : 


OM V 

Kyiy±UU’> 


OM"  OM' 


On  pose  : 

A = n.  OM'uu , M = n.OM'uv,  N = n.OM'vv 

9.  Calculer  les  coefficients  L , M et  N. 

10.  En  déduire  l’expression  des  courbures  moyenne  et  totale. 


Problème  n° 7 : Soit  la  surface  (r)  définie  paramétriquement  par  : 


M (u,  v ) < 


X = thu.cosv 
Y = thu.sinv 

1 U 

Z=7h~u  + LOSth2 


avec  chu  et  thu  sont  respectivement  le  cosinus  et  la  tangente  hyperboliques  définies  par  : 


U | „—u 


chu  = 


e + e 


U | _ — U 


thu  = 


e + e 


U SD  — U 


e“  — e 


1.  Donner  le  domaine  de  définition  de  la  surface  (r). 

2.  Calculer  les  composantes  des  vecteurs  OMu  et  OMv. 
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3.  Calculer  les  coefficients  E,  F,  G de  la  première  forme  fondamentale  de  la  surface 


4.  En  déduire  l’expression  de  ds2. 

5.  Les  coordonnées  (u,v)  sont-elles  orthogonales?  symétriques? 

6.  Calculer  un  vecteur  unitaire  normal  n de  (r). 

7.  Calculer  les  vecteurs  : 


OM"au,  ONT' 


OM" 

vv 


On  pose  : 

L = n.OM'uu,  M = n.OM'uv , N = n.OM'vv 

8.  Calculer  les  coefficients  L,  M et  N. 

9.  Déterminer  les  coubures  moyenne  et  totale. 


Problème  n° 8 : Montrer  que  les  courbures  totale  K et  moyenne  H en  un  point 
M(x,y,z ) d’une  surface  paramétrée  par  z = f(x,y),  où  / est  une  fonction  lisse,  sont 
données  par  : 

fil  fil  _ fl  12 

jy-  _ JxxJyy  Jxy 

~ (1  +f'X2+f!y2)2 
et  : 


H = 


a + f'x2)f'L  - v'j'r'  + a + s^r: 


yy 


(i  + f'x2  + K2 


Problème  n° 8 : Soit  (E)  une  surface  de  M3  paramétrée  par  OM(u,v ) telle  que  sa 
première  forme  fondamentale  s’écrit  : ds2  = Edu 2 + 2 Fdudv  + Gdv2 

1.  Montrer  que  les  conditions  suivantes  sont  équivalentes  : 

(Æ  _ dG  _ 

dv  du  1 

d2OM 

ii)  - Le  vecteur  — est  parallèle  au  vecteur  normal  N à la  surface, 

auov 

iii)  - Les  côtés  opposés  de  tout  quadrilatère  curviligne  formés  par  les  courbes  co- 
ordonnées (u,  v ) ont  même  longueurs. 

2.  Quand  ces  conditions  sont  satisfaites,  on  dit  que  les  courbes  coordonnées  de  (E) 
forment  un  réseau  de  Tchebychev.^  Montrer  que  dans  ce  cas,  on  peut  paramétrer  la 
surface  par  (ü,  v)  telle  que  ds2  s’écrit  : 

ds2  = dû2  + 2 cosOdüdv  + dv2 


3.  Pafnouti  Tchebychev  (1821  - 1894  ) : Mathématicien  russe. 
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où  6 est  une  fonction  de  (ü,v).  Montrer  que  9 est  l’angle  entre  les  courbes  coordonnées 
ü,  v. 

3.  Montrer  que  l’expression  de  la  courbure  totale  est  donnée  par  : 


K = 


i d2e 

sinO  dûdv 


4.  On  pose  : 


û = û + v 


V = u — V 

Montrer  que  ds2  s’écrit  avec  les  nouvelles  variables  ( û , v ) : 

O 0 0 0 0 

ds  = cos  u dû  + sin  lj  dû 
avec  lj  = 9/2.  (A. N.  Pressley,  2010) 
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3.4  La  Géométrie  de  l’Ellipse  et  de  l’Ellipsoïde 

Exercice  n°l  : A partir  de  la  définition  géométrique  de  l’ellipse  donnée  par  : 

MF  + MF'  = constante  = 2 a 
retrouver  l’expression  de  l’équation  cartésienne  de  l’ellipse. 


Exercice  n° 2 : Montrer  la  formule  très  utilisée  en  géodésie  : 

d(Ncosp) 


dp 


= —psmp 


avec  N et  p les  deux  rayons  de  courbures  principaux  de  l’ellipsoide  de  révolution  donnés 
respectivement  par  : 

a, 

N = 


et  : 


p = N = 


y/l  — e2sin2ip 


a(l  — e2) 


(1  — e2sin2p)yj  1 — e2sin2p 


Problème  n°l  : A partir  des  équations  de  l’ellipsoide  de  révolution  : 


M = 


X = N cospcosX 
Y = NcospsinX 
Z = N(  1 — e2)sinp 


1.  Calculer  les  vecteurs  : 


2.  Calculer  les  coefficients  : 


dM  d M 

~dX'~d p 


dM  dM  _ dM  dM  dM  dM 

~~dx'~dx,  ~ dx  ' dp  ' ^~lk/'~d^ 

Démontrer  que  l’expression  de  la  première  forme  fondamentale  s’écrit  : 


ds 2 = p2dip2  + N2cos2(pdX 2 


3.  Calculer  le  vecteur  normal  n : 


dM  dM 

dX  ^ dtp 


1 

~dM  dM 
dX  ^ dp 


n = 
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4.  Calculer  les  vecteurs  : 


5.  Déterminer  les  coefficients  : 


d2M  d2M  d2M 


dX2  ' dXdip'  dip2 


d2M  d2M  d2M 

L = n.  ^ , M = n.  n ^ , N = n.- 


d\2  ’ ' dX dp> 

6.  Ecrire  la  deuxième  forme  fondamentale  <f>(A,(/?). 

7.  En  appliquant  la  formule  du  cours,  Montrer  que  : 

N(<p)  = 


d2<p 


i/l  — e2sin2p> 


est  le  rayon  de  courbure  de  la  section  normale  au  point  M perpendiculaire  au  plan  de  la 
méridienne  de  l’ellipsoide  de  révolution. 

8.  En  posant  : 


dL  = 


pdp> 

Ncosip 


En  déduire  que  ds 2 s’écrit  : 


ds 2 = N~cos2(p(d£~  + dX2) 


9.  Montrer  que  C est  donnée  par  : 

£(»  = Log 


7T  (fi 

4 + 2 


e / 1 + esirup 

~ 3L°9  ï : — 

2 \ 1 — esirup 


Problème  n° 2 : Sur  l’ellipsoide,  on  note  ip  la  latitude  géodésique  et  ip  la  latitude 
réduite. 

1.  Calculer  p le  rayon  de  courbure  de  l’ellipse  méridienne  en  fonction  de  'î/c 

2.  Exprimer  l’aplatissement  de  l’ellipsoide  en  fonction  des  valeurs  de  p au  pôle  et  à 
l’équateur. 

3.  On  mesure  la  longueur  d’un  arc  de  méridien  d’un  degré  à la  fois  au  pôle  et  à l’équa- 
teur. On  trouve  respectivement  111695  m et  110  573  m.  En  déduire  l’aplatissement. 

Problème  n°3  : On  donne  les  coordonnées  tridimensionnelles  suivantes  d’un  point 
M : 

M = (X,  Y,  Z)  = (4  300  244.860  m,  1 062  094.681  m,  4 574  775.629  m) 

Les  paramètres  de  l’ellipsoide  de  référence  sont  a = 6 378  137.00  m,  e2  = 0.006  694  38. 
1.  Calculer  le  demi-petit  axe  b. 
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2.  Calculer  l’aplatissement. 

3.  Calculer  les  coordonnées  géodésiques  (</?,  A,  he ) du  point  M.  ip  e t À seront  calculées 
en  grades  avec  cinq  chiffres  après  la  virgule. 

Problème  n°4  : Soit  S (a,  e ) un  ellipsoide  de  révolution  où  a,  e sont  respectivement  le 
demi-grand  axe  et  la  première  excentricité.  ( g ) une  géodésique  partant  d’un  point  E(tp  = 
0,  A^;)  sur  l’équateur  et  d’azimut  Aze • A cette  géodésique,  on  lui  fait  correspondre  une 
géodésique  ( g ')  sur  la  sphère  S 2 dite  de  Jacobi|^]de  rayon  a,  ayant  le  même  azimut  Aze 
au  point  E\gJ  = 0,  A^;).  De  même  au  point  M(ip,  A)  de  la  géodésique  ( g ) de  l’ellipsoide, 
on  lui  fait  correspondre  le  point  M\gJ , A;)  de  (g1)  de  S 2 tel  qu’il  y a conservation  des 
azimuts. 


P 


Figure  3.4.1  - La  Correspondance  de  la  sphère  de  Jacobi 


1.  Ecrire  l’équation  de  Clairaut  pour  la  géodésique  (g). 

2.  On  note  r'  le  rayon  du  parallèle  passant  par  M'  de  la  géodésique  {g').  Ecrire  de 
même  l’équation  de  Clairaut  pour  la  géodésique  ( g 

3.  Montrer  que  p et  gJ  vérifient  : 

Ncosp  = acosgJ 


et  en  déduire  que  gJ  est  la  latitude  paramétrique  de  AI. 

4.  Ecrire  les  expressions  de  tgAzg  et  tgAzg’  respectivement  sur  (g)  et  {g'). 

5.  Montrer  que  : 


dX 


pdg) 

adg>' 


d\' 


En  déduire  que  : 


d\  = a/1  — e2cos2>g>'d\' 


6.  En  intégrant  l’équation  précédente,  montrer  qu’on  obtient  : 


A — A e — 


e2cos2g>'d\' 


4.  Cari  Gustav  Jacob  Jacobi  (1804-1851)  : Mathématicien  allemand. 
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avec  À > \e  et  A'  est  comptée  à partir  de  A E- 

7.  En  écrivant  yQ  — e2cos2(p'  = 1 — ^ cos 2ip'  + o(e4)  où  o(e4)  est  un  infiniment  petit 
d’ordre  4 en  e dont  on  néglige,  écrire  l’intégrale  précédente  entre  XE  et  A^  + A. 

8.  Comme  ( g ')  est  une  géodésique  de  la  sphère,  on  démontre  que  : 


cos2^'dX! 


siuAze  , , 

as 

a 


où  ds'  est  l’élément  différentiel  de  l’abscisse  curviligne  sur  la  géodésique  (un  grand  cercle). 
Alors  en  posant  s'  = 0 au  point  E' , montrer  que  l’équation  précédente  s’écrit  sous  la 
forme  : 


A = XE  + A'  - 


e2sinAz. 


2 a 


e I ds' 
o 


9.  On  suppose  que  la  géodésique  ( g ')  coupe  une  première  fois  le  plan  de  l’équateur  en  un 
point  F' , montrer  qu’on  obtient  : 


A e = tt 
s'  = tt  a 

, , e27 TsinAzE 

\f  = \e  + tt 

10.  La  géodésique  (g1)  partant  de  F'  a pour  azimut  n — AzE,  elle  coupe  une  deuxième 
fois  l’équateur  au  point  E1 , mais  la  géodésique  (g)  sur  l’ellipsoide  coupe  une  deuxième 
fois  le  plan  de  l’équateur  au  point  correspondant  à H dont  la  longitude  est  A h-  Montrer 
que  A h est  donnée  par  : 

, , e27 TsinAzE  e2Trsin(TT  — Aze)  „ 9 

\H  = \e  + 2tt FL  = \E  + 2tt  - e~TrsinAzE 

Quelle  conclusion  a-t-on  sur  les  lignes  géodésiques  de  l’ellipsoide  de  révolution. 

Problème  n° 5 : Un  point  M de  la  surface  d’une  sphère  (S)  de  rayon  R,  a pour 
coordonnées  (A,  Y,  Z)  dans  un  repère  orthonormé  : 

M = (A,  Y,  Z)  = (Rcosp.cos A,  Rcoscp.sinX , Rsirup) 

1.  Montrer  qu’un  vecteur  normal  unitaire  n à ( S ) en  M est  : 

n = ( cosip.cosX , cosip.sinX,  simp)T 

2.  Soit  (C)  le  grand  cercle  passant  par  le  point  A(R,  0,0)  et  d’azimut  AzE-  Le  point 
M peut  être  décrit  par  son  abscisse  curviligne  s mesurant  l’arc  AM.  On  note  par  eu 
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représente  l’angle  au  centre  de  l’arc  AM.  Utilisant  la  trigonométrie  sphérique,  montrer 
que  : 

costp.sinX  = sinoj.sinAzE 

3.  En  utilisant  la  formule  fondamentale  de  la  trigonométrie  sphérique  dans  le  triangle 
APM,  montrer  qu’on  a les  deux  relations  : 


cosoj  = costp.cosX 
sirup  = sinuj.cosAzE 


4.  En  déduire  que  les  coordonnées  de  M s’écrivent  en  fonction  de  s comme  suit  : 


M < 


X = R.cos(s/R) 

Y = RsinAzEsin(s  / R) 
Z = RcosAzEsin(s  / R ) 


5.  Calculer  les  vecteurs  T et  N du  repère  de  Frenêt.  En  déduire  les  composantes  de  N 
en  fonction  de  oj. 

6.  Montrer  que  les  vecteurs  N et  n sont  parallèles. 

7.  Justifier  que  les  géodésiques  de  la  sphère  sont  les  grands  cercles. 


Problème  n°6  : Soit  le  tore  T défini  par  les  équations  suivantes  : 


M(tp,X) 


x = (a  + Rcostp)cosX 
< y = (a  + Rcostp)sinX 
z = Rsinip 


où  a,  R deux  constantes  positives  avec  a > R,  (tp,  A)  G [0,  2n\  X [0,  2tt\. 

1.  Calculer  la  première  forme  fondamentale  ds2. 

2.  Avec  les  notations  usuelles,  on  pose  : 


dE  . dE  , dF  . 

= E = E[  — = F 

<P'  d\  A’  P),n  v 


dtp 


dtp 

dG 


dF  _ , dG  _ , _ , 

~dX  ~ ÔX~Gx 

Utilisant  les  équations  des  géodésiques  du  cours,  montrer  que  les  équations  des  géodé- 
siques du  tore  sont  : 

dtp  dX  o d2  X 

— 2RsiTitp[(i  + Rcostp)— — - — |-  (a  + Rcostp)  = 0 
ds  ds  dsz 


Rsirup(a  + Rcostp)  ( — ) + R2-j-^  = 0 


>d2tp 
ds 2 
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3.  Montrer  que  la  première  équation  ci-dessus  donne  : 

(a  + Rcosip)2  — = C = cte 
as 


Montrer  qu’on  retrouve  l’équation  de  Clairaut  avec  C = (a  + R)sinAze  où  Aze  est 
l’azimut  de  départ  au  point  Mo(<p  = 0,  Ao). 

4.  On  suppose  au  point  Mo,  la  géodésique  a pour  azimut  Aze  tel  que  : 


0 < Aze  < 


TT 

2 


Montrer  que  la  deuxième  équation  des  géodésiques  s’écrit  en  utilisant  le  résultat  précé- 
dent : 

d2ip  C 2 sirup 
ds2  R (a  + Rcostp )3 

5.  Montrer  qu’on  arrive  à : 


dtp 

ds 


Ml  =i_ 


C 2 


R2  (a  + Rcostp )2 


> 0 


où  l est  une  constante  d’intégration. 
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3.5  Les  Systèmes  Géodésiques 

Exercice  n°  1 : Donner  l’expression  des  composantes  du  gradient  en  coordonnées 
cylindriques. 

Exercice  n° 2 : On  donne  l’expression  scalaire  d’une  fonction  V(x,y,z)  par  : 

T//  \ CLX2  + y2  1 2 , 2 2. 

V(x,y,z)  = 2 + (x  + 2/  ) 

2 

1.  Calculer  les  composantes  du  vecteur  gradV  dans  un  domaine  de  R3  où  z ^ 0. 

Problème  n°  1 : Soit  un  point  A(<p,  À)  sur  un  ellipsoide  de  révolution  associé  à un  ré- 
férentiel géocentrique  donné  7Z.  On  considère  le  repère  orthonormé  local  en  A ( e\ , e^,  en) 
défini  dans  la  base  orthonormée  (i.  j,  k ) de  1Z  où  e\  est  tangent  au  parallèle  passant  par 
A et  dirigé  vers  l’Est,  tangent  à la  méridienne,  dirigé  vers  le  nord  et  en  porté  par  la 
normale  à l’ ellipsoide  dirigé  vers  le  zénith. 

1.  Exprimer  les  vecteurs  de  la  base  (e,\,  e^,  en)  dans  la  base  ( i,j,k ) de  7 Z. 

2.  Exprimer  les  vecteurs  i , j et  k dans  la  base  (e\,e,p,en). 

3.  Calculer  de\,dev  et  den  dans  la  base  (■ i,j,k ). 

4.  En  adoptant  une  écriture  matricielle,  montrer  que  : 


( de\\ 

( 0 

sirupdX 

—cos(pdX\ 

de^ 

= 

—sirupdX 

0 

—dtp 

\dCn  j 

\ cos(fd\ 

dtp 

o ) 

\en) 

Problème  n° 2 : On  définit  dans  R3  un  point  M par  ses  coordonnées  ellipsoidiques 
de  Jacobi  (cfr,  À,  u)  comme  suit  : 


M < 


x = Vu2  + e2.cos(f>cos\ 
y = \Ju2  + e2  .cos(f)sin\ 
z = u.siruj) 


avec  : e2  = \J a2  — b2 , (j)  G [— 7r/2,  7t/2],  A G [0,  27t]  et  u G]0,  +oo[,  a,  b deux  constantes 
réelles  telles  que  a > b > 0. 

1.  Montrer  que  le  point  M appartient  à un  ellipsoide  de  révolution  en  précisant  ses 


demi-axes. 
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2.  Calculer  ds 2 et  montrer  qu’il  s’écrit  sous  la  forme  : 


ds 2 = ( d(j),d\,du).G . 


( d<f>\ 

dX 

\du  ) 


avec  G donnée  par  : 


G = (gij)  = 


( u2  + e2sin24>  0 

0 ( u 2 + e2)co.s2(f> 

\ 0 0 


0 \ 

0 

u 2 + e2sin2c 


u2  + e2  / 

3.  Sachant  que  l’expression  du  laplacien  d’une  fonction  scalaire  V en  coordonnées  de 
Jacobi  est  exprimée  par  : 


AT,  = ±{JLf>@W\  dfy/gâV\  d_  f y/g 

y/g\d4>\gu' d(t>  ) d\\g22  d\  J du  \g33' du  J j 

où  g est  le  déterminant  de  la  matrice  G,  donner  l’expression  de  AV. 

4.  Calculer  AV  sachant  que  V est  donnée  par  : 

jr  f,  \ GM  ,e.  ^ 2 2e!  fi  ^ 2 A , 1 2/  2 . 2\  2 i 

V (<p,  u)  = Arctg — I — au  — 1 cos  cp  + -u  (u  + e jcos  ç> 

e u 3 qo  \ 2 J 2 

avec  G,  M et  u des  constantes  et  : 


q = q(u ) = 


1 


1 + 3- 


u 


u 


Arctg 3 


u 


Qo  = q(u  = b)  = - 


b2 


1 + 3-^  Arctg-  — 3- 
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3.6  Les  Réductions  des  Distances 

Exercice  n°  1 : On  a mesuré  une  distance  suivant  la  pente  Dp  = 20130.858  m entre 
deux  points  A et  B avec  Ha  = 235.07m,  Hp  = 507.75  m,  on  prendra  comme  rayon 
terrestre  R = 6378  km .. 

1.  Calculer  la  distance  suivant  l’ellipsoide  : 

- en  utilisant  les  différentes  corrections, 

- en  utilisant  la  formule  rigoureuse. 

2.  En  prenant  la  valeur  de  la  formule  rigoureuse  et  sachant  que  le  module  linéaire  m 
vaut  0.999  850  371,  calculer  la  distance  réduite  au  plan  de  la  représentation  plane  utilisée. 

Exercice  n°2  : Entre  2 points  A ( Ha  = 128.26  m ) et  B ( Hp  = 231.84  m),  la 
distance  Dp  suivant  la  pente  est  égale  à 15  498.823  m.  Soit  Dq  la  distance  corde  au 
niveau  de  la  surface  de  référence.  L’angle  de  site  observé  en  A en  direction  de  B est 
i = 0.3523  gr. 

1.  Calculer  la  valeur  de  Dq  en  utilisant  la  formule  rigoureuse. 

2.  Calculer  Dq  par  les  corrections. 

3.  En  adoptant  la  moyenne  des  deux  méthodes,  calculer  la  distance  De  réduite  à la 
surface  de  référence. 

4.  Le  module  linéaire  de  la  représentation  plane  Lambert  Sud  utilisée  est  de  0.999  648  744. 
calculer  alors  la  distance  Dr  réduite  au  plan  de  la  représentation. 

Exercice  n°3  : On  a mesuré  une  distance  suivant  la  pente  entre  les  points  A (Ha  = 
1319.79  m)  et  B (Hp  = 1025.34  m)  avec  Dp  = 16  483.873  m. 

1.  Calculer  la  distance  De  distance  réduite  à l’ellipsoide  de  référence  par  la  formule 
rigoureuse,  on  prendra  le  rayon  de  la  Terre  R = 6378  km. 

2.  Calculer  la  distance  Dr  réduite  à la  représentation  plane  Lambert  si  l’altération 
linéaire  de  la  zone  est  de  —14  cm/ km. 
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3.7  Les  Représentations  Planes 


Problème  n°  1 : Soit  S2  la  sphère  de  rayon  R , au  point  P(<p,  A)  on  lui  fait  corres- 
pondre le  point  p(X,  Y)  du  plan  OXY  par  la  représentation  plane  suivante  définie  par 
les  formules  : 


p(X,Y)  = 


X = 2R.tg{j  - ^ ).sin\ 
Y = -2 R.tg(j  - |).co.sA 


1.  Montrer  que  l’image  d’un  méridien  (À  — constante  ) est  une  droite  dont  on  donne 
l’équation. 

2.  Montrer  que  l’image  d’un  parallèle  (<p  — constante  ) est  un  cercle  dont  on  précise 
l’équation. 

3.  En  utilisant  le  lennne  de  Tissot,  déterminer  les  directions  principales. 

4.  Soit  dS  la  longueur  infinitésimale  correspondante  sur  le  plan,  calculer  dS. 

5.  Sachant  que  sur  la  sphère  ds2  = R2dp 2 + R2 cos2 p .d\2 , calculer  le  module  linéaire 


m. 

6.  En  déduire  le  module  linéaire  rri\  le  long  du  méridien. 

7.  En  déduire  le  module  linéaire  m?  le  long  d’un  parallèle. 

8.  Comparer  m\  et  rri2.  Conclure  sur  la  conformité  ou  la  non  conformité  de  la  repré- 
sentation plane. 


Problème  n° 2 : Soit  S la  sphère  de  rayon  R,  au  point  P(p,  A)  on  lui  fait  corres- 
pondre le  point  p(X,  Y)  du  plan  OXY  par  la  représentation  plane  suivante  définie  par 
les  formules  : 


p(X,Y) 


X = R.  A 

Y = R.Logtg(j  + |) 


où  Log  désigne  le  logarithme  népérien. 

1.  Quelles  sont  les  images  des  méridiens  (A  = constante)  et  des  parallèles  (<p  = 
constante). 

2.  Soit  dS  la  longueur  infinitésimale  correspondante  sur  le  plan,  calculer  dS  en  fonc- 
tion de  p et  de  A et  calculer  le  module  linéaire  m. 

3.  En  déduire  les  modules  linéaires  m\  le  long  du  méridien  et  m2  le  long  du  parallèle. 

4.  Comparer  mi  et  777.2  et  conclure  sur  la  conformité  ou  la  non  conformité  de  la 
représentation  plane. 

5.  On  suppose  que  P décrit  sur  la  surface  S une  courbe  (7)  telle  que  (p  et  A sont 
liées  par  la  relation  : tgip  = sinX.  Pour  <p  = 0 gr,  2 gr , 4 gr,  6 gr,  8 gr  et  10  gr,  dresser  un 
tableau  donnant  les  valeurs  de  A correspondantes. 


44 


CHAPITRE  3.  EXERCICES  ET  PROBLEMES 


6.  Sachant  que  R = 1000  m,  calculer  les  coordonnées  (X,Y)  de  la  représentation 
plane  donnée  ci-dessus  pour  les  valeurs  de  <p  et  A de  la  question  5. 

7.  Rapporter  à l’échelle  1/100  sur  le  plan  OXY,  les  positions  (X,Y)  des  points.  Que 
pensez- vous  de  l’image  de  la  courbe  (7). 

Problème  n° 3 : Sur  une  sphère  de  rayon  unité,  modèle  de  la  terre,  on  désigne  : 

- par  p le  pôle  nord, 

- par  (C)  un  grand  cercle  qui  coupe  l’équateur  au  point  i de  longitude  nulle, 

- par  q le  pôle  de  ce  grand  cercle,  de  latitude  ipQ  positive, 

- par  uj  et  h respectivement  les  points  d’intersection  de  ( C ) et  du  méridien  de  q et 
du  grand  cercle  issu  de  g,  passant  par  le  point  a(ip,  A). 

On  pose  : ujh  = x,  ha  = y 

1.  q est  le  pivot  d’une  représentation  cylindrique  conforme  oblique  tangente,  dont 
(C)  est  le  ”pseudo-équateur”.  Le  plan  est  rapporté  aux  axes  DAT,  LlY  images  respectives 
de  (C)  et  du  grand  cercle  upq.  Exprimer  en  fonction  de  (/?,  A et  içq  les  coordonnées  X,  Y 
du  point  A image  de  a (vérifier  que  pour  ipQ  = 0,  on  retrouve  les  expressions  de  X,  Y 
d’une  représentation  transverse). 

2.  Montrer  que  l’équation  de  l’image  plane  du  parallèle  de  latitude  ip$  peut  s’écrire  : 

e1  cosX  = tgtpo 

Indications  : b désignant  un  point  de  latitude  ipo,  le  triangle  pqb  est  isocèle,  décom- 
poser ce  triangle  en  deux  triangles  rectangles  égaux.  Etudier  qualitativement  les  images 
des  autres  parallèles. 

3.  Montrer  que  l’image  plan  de  l’équateur  a pour  équation  : 

cosX  + tgcpo-shY  = 0 

Ecrire  d’une  manière  analogue,  l’équation  de  l’image  du  méridien  A = 0. 

4.  Exprimer  le  gisement  du  méridien  en  fonction  de  ip,  A et  ipç,.  Déterminer  la  valeur 
du  module  linéaire,  en  particulier  en  p,  en  un  point  de  l’équateur,  en  un  point  du  méri- 
dien origine. 

Problème  n°4  : Etude  de  la  représentation  conforme  d’une  sphère  de  rayon  unité 
dite  représentation  de  Littrow]^]  définie  par  : 

Z = sinz 


5.  En  hommage  à Joseph  Johann  Littrow  (1781-1840)  astronome  autrichien. 
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avec  z = X + iLetZ  = X + ïY . 

1.  Préciser  le  canevas,  les  images  des  méridiens  et  celle  de  l’équateur. 

2.  Vérifier  que  les  points  / et  f (<p  = 0,  A = ±vr/2)  sont  des  points  singuliers. 

3.  Etudier  les  images  plans  des  cercles  de  diamètre  f f et  des  petits  cercles  orthogo- 
naux. 

4.  Soit  s le  point  (p  = ipo,  A = 0).  On  appelle  segment  capable  sphérique  l’ensemble 
des  points  b tels  que  l’angle  bp , bs  = a.  Quelle  est  l’image  plane  de  cette  courbe  dans 
cette  représentation  plane. 

Problème  n° 5 : Soit  l’application  F(u,v ) : R2  — > R3\(0,0, 1)  définie  par  : 

2 u 

2 1 2 i T" 

uz  + vz  + 1 

2v 

OM{u,v)  = F(u,v)  y=  u2  + v2  + 1 

u2  + v2  — 1 

2 = TT 

UZ  + V2,  + 1 

1.  Calculer  la  forme  fondamentale  ds2 . 

2.  Montrer  que  OM(u,v)  appartient  à la  sphère  S2  d’équation  x2  + y2  + z2  = 1. 

3.  Calculer  u,  v en  fonction  de  x.  y et  z. 

4.  Soit  le  point  N (0,  0, 1)  de  S2,  calculer  les  coordonnées  ( X , Y)  du  point  p intersection 
de  la  droite  NM  avec  la  plan  z = 0 en  fonction  de  x,y  et  z. 

5.  Soit  a l’application  R3\(0,  0, 1)  — > R2  : (x,  y,  z)  — > ( X , Y)  = (X(x,  y , z),  Y (x,  y,  z )). 
Montrer  que  (a  o F)(u , v ) = a(F(u,  v ))  = ( u , v).  En  déduire  que  F = cr'1. 

6.  Trouver  le  rapport  de  ce  problème  avec  le  problème  ??. 

Problème  n°6  : Soit  un  ellipsoide  de  révolution  E(a,e)  avec  a et  e respectivement 
le  demi-grand  axe  de  l’ellipsoide  de  révolution  et  e la  première  excentricité.  Soit  S2  une 
sphère  de  rayon  R.  On  veut  étudier  le  passage  suivant  : 

p((p,  À)  de  l’ellipsoide  E P(ip,  A)  de  la  sphère  S2 

1.  Exprimer  m le  module  linéaire  de  cette  représentation. 

2.  On  pose  : 

z = jC.  -{-  zA,  Z — L -\-  iA 

C est  la  latitude  isométrique  de  l’ellipsoide  de  révolution  et  L la  latitude  de  Mercator. 
Une  transformation  conforme  entre  E et  S2  est  donnée  par  Z = f(z)  où  / est  une 


CHAPITRE  3.  EXERCICES  ET  PROBLEMES 


fonction  analytique.  On  propose  le  cas  le  plus  simple  à savoir  : 


Z = az  + /3 
et  = Cr  ic2 
P = b\  + ib2 


les  ci,  C2,  b\ , b2  sont  des  constantes  réelles.  Donner  les  expressions  de  L et  A en  fonction 
de  C et  À. 

3.  On  veut  que  repésentation  transforme  les  méridiens  et  les  parallèles  de  l’ellipsoide 
respectivement  en  méridiens  et  parallèles  de  la  sphère  et  que  l’image  du  méridien  origine 
A = 0 soit  le  méridien  origine  de  la  sphère  A = 0.  Montrer  que  c2  = b2  = 0 et  L = 
c\L  + b\,  A = ciÀ. 

4.  Pour  avoir  la  même  orientation  en  longitude,  on  prendra  ci  > 0.  On  cherchera  la 
transformation  à déformation  minimale  autour  d’un  parallèle  ip  = <pq  tel  que  le  parallèle 

ip  = tpo  est  automécoique  et  le  module  linéaire  m est  stationnaire  pour  (p  = (po,  c’est-à- 

(din  \ 

) = 0,  en  plus  on  considère  aussi  la  condition  : 

d(f  J 


Pour  faciliter  les  notations,  on  prendra  b = b\,  c = c±.  Montrer  que  la  relation  liant  (po 
et  ’I'o  est  : 


tg'P  o = tg<po\ 


1 — e2sin2ipo 


5.  Déterminer  les  constantes  6,  c et  R en  fonction  de  ip o et  'Pq  telles  que  les  conditions 
ci-dessus  soient  vérifiées. 

6.  Montrer  que  l’expression  du  développement  limité  de  m(<p)  de  part  et  d’autre  du 
parallèle  <pç>  est  donnée  par  : 

2e2  (1  — e2)sirupocosipo  3 4 

m(v)  = 1 - 3(l-A,mW  (y  ~ W)  + °i(V  ~ W)  ] 

7.  On  fait  intervenir  la  deuxième  excentricité  e' , Montrer  que  m(ip)  s’écrit  : 

2e,2sirup0cosip0  3 4 

m(y>  = 1~3(l  + e-Wyo)4*’-^  +o((v-vo)) 


Problème  n° 7 : Soit  £(a,b ) un  ellipsoide  de  référence  de  paramètres  a et  e respec- 
tivement le  demi-grand  axe  et  la  première  excentricité.  On  considère  une  représentation 
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plane  V de  £ vers  le  plan  (O,  X,  Y).  On  pose  : 

z = A + iC 
Z = X + iY  = Z(z) 

avec  C la  latitude  isoparamétrique. 

1.  Ecrire  les  expressions  du  carré  des  éléments  infinitésimaux  de  longueur  sur  l'ellip- 
soïde et  le  plan.  En  déduire  le  module  linéaire  m. 

dZ 

2.  On  pose  C = — — . Si  7 est  le  gisement  de  l’image  du  méridien  passant  par  le  point 

oz 

7 r 

z = (À,  £),  montrer  que  arg(Q  = — — 7. 

3.  On  cherche  une  représentation  plane  du  type  Z = a + {3z  + wz2  où  a,  fi  et  zu  des 
constantes  complexes.  On  impose  les  conditions  suivantes  : 

- pour  z = 0,  Z = 0, 

- l’axe  des  Y coincide  avec  le  méridien  à l’origine. 

Montrer  que  !Ze(f3)  = 0. 

4.  En  déduire  que  Z s’écrit  : 

Z = i/3\z  + (vj  1 + ivo2)z2 


avec  ^1,071,^372  sont  des  réels. 
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CHAPITRE  3.  EXERCICES  ET  PROBLEMES 


3.8  La  Représentation  Lambert 

Exercice  n°l  : En  un  point  A de  coordonnées  géodésiques  <p  = 40.9193  gr  et  À = 
11.96565?’  à l’Est  de  Greenwich,  on  vise  un  point  B. 

1.  Dans  quelle  zone  de  Lambert  Tunisie  se  trouve  le  point  A ? Calculer  ses  coordonnées 
planes  (X,Y). 

2.  L’azimut  géodésique  de  la  direction  AB  est  Azg  = 55.7631  gr.  Sachant  que  Dv  = 
1.52  dm.gr , calculer  G le  gisement  de  la  direction  AB. 

3.  La  distance  AB  réduite  à l’ellipsoide  de  référence  est  De  = 5421.32  m..  Sachant 
que  l’altération  linéaire  dans  la  région  des  points  A et  R vaut  —9  cm/ km,  calculer  la 
distance  AB  réduite  au  plan. 

Exercice  n° 2 : D’après  les  coordonnées  de  deux  points  A et  R vous  trouvez  la 
distance  AB  = 5427.380  m.  Sachant  que  : 

a - l’altération  linéaire  de  la  représentation  dans  la  région  de  AB  vaut  +8.10“5, 

b - les  altitudes  des  points  A et  R sont  : Ha  = 1000.00???.  et  Hb  = 1200.00  m.  Cal- 
culer la  distance  suivant  la  pente  Dp  entre  les  points  A et  R matérialisés  sur  le  terrain. 

Problème  n°l  : On  a mesuré  une  distance  suivant  la  pente  Dp  = 20130.858???. 
entre  deux  points  A et  R avec  Ha  = 235.07???,  Hp  = 507.75  m,  on  prendra  comme 
rayon  terrestre  R = 6378  km. 

1.  Calculer  la  distance  De  suivant  l’ellipsoide  en  utilisant  la  formule  rigoureuse. 

2.  Sachant  que  le  module  linéaire  m vaut  0.999  850  371,  calculer  la  distance  Dr  réduite 
au  plan  de  la  représentation  plane  utilisée. 

3.  Les  coordonnées  géodésiques  du  point  A sont  : <p  = 10.7245  3 gr,  À = 41.4490  3 gr. 
Par  des  observations  astronomiques,  on  a déterminé  les  coordonnées  astronomiques  (pa  = 
10.7257  4 gr  et  Aa  = 41.4505  2 gr  du  point  A et  l’azimut  astronomique  de  la  direction 
AB  soit  Aza  = 89.6849  9 gr.  Transformer  l’azimut  astronomique  de  la  direction  AB  en 
azimut  géodésique  en  utilisant  l’équation  de  Laplace  donnée  par  : 

Azg  = Aza  + (A  — A a).sirup 

4.  Calculer  le  gisement  G de  la  direction  AB  sachant  que  Dv  = 0.0018  8 gr. 

5.  Les  coordonnées  Lambert  Nord  Tunisie  de  A sont  X = 478  022.43???.  et  Y = 
444  702.22???.  Déterminer  alors  les  coordonnées  de  R. 

6.  Calculer  l’azimut  de  R vers  A sachant  qu’on  néglige  la  correction  de  la  corde  de 
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la  direction  B A et  que  Xp  = 10.9288  4 gr. 

Problème  n° 2 : On  a mesuré  une  distance  suivant  la  pente  entre  les  points  A(Ha  = 
1319.79  m)  et  B(Hb  = 1025.34m)  avec  Dp  = 16483.873  m. 

1.  Calculer  la  distance  De  distance  réduite  à l’ellipsoide  de  référence  par  la  formule 
rigoureuse,  on  prendra  le  rayon  de  la  Terre  R = 6378  km. 

2.  Calculer  la  distance  Dr  réduite  à la  représentation  plane  Lambert  si  l’altération 
linéaire  de  la  zone  est  de  —14  cm/km. 

3.  La  direction  AB  a un  azimut  géodésique  Azg  = 297.5622  5 gr.  Donner  l’expression 
du  gisement  G de  AB  en  fonction  de  Azg, 'y  la  convergence  des  méridiens  et  Du  la 
correction  de  la  corde,  sachant  que  la  représentation  plane  utilisée  est  le  Lambert  Sud 
Tunisie  et  que  le  point  A est  au  nord  du  parallèle  origine. 

4.  On  donne  Dv  = —13.7  dmgr  et  À = 9.3474  734  gr  la  longitude  de  A,  calculer  G. 

5.  En  déduire  les  coordonnées  (Xp,  Yp)  de  B si  Xa  = 363  044.79  m et  Y a = 407  020.09  m. 

6.  Déterminer  les  coordonnées  géographiques  {ip,  A)  de  B. 

On  rappelle  que  : a = 6 378  249.20  m.  et  e2  = 0.006  803  487  7. 
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CHAPITRE  3.  EXERCICES  ET  PROBLEMES 


3.9  La  Représentation  UTM 

Exercice  n°  1 : Dans  cet  exercice,  on  voudrait  justifier  l’arrêt  à l’ordre  8 de  l’expres- 
sion de  Y (UTM)  en  fonction  de  À.  On  donne  : p = 40.00  gr  et  a = 6 378  249.20  m,  e1 2  = 

0.006  803  4877. 

1.  Calculer  numériquement  e/2,r/2,f2 * *  = tg2 *p  et  N(p). 

2.  Calculer  numériquement  le  coefficient  as  de  (??). 

3.  On  donne  À = 1.235  46  gr,  calculer  as  A8  et  conclure. 


Problème  n°l  : Soit  le  point  A de  coordonnées  géodésiques  : p = 40.9193 gr  et 
À = 11.9656  gr  à l’Est  de  Greenwich.  On  considère  la  représentation  plane  UTM  tronquée 
suivante,  de  méridien  central  Ào  = 9°  définie  par  les  formules  : 

X = ai.(À  — Ào)  + 03. (A  — Ao)3 
Y = g(p)  + a 2.(A- A0)2 
où  p,  À et  Ao  sont  exprimées  en  rd,  avec  : 

ai  = N(p).cosp 
a\  . 

ü2  = —.smp 


03  = 


a\cos  p 
6 

N(p)  = 


(1  — tg2p  + e2  .cos2p) 
a 


eA  .sin- p 


g(p)  = a(l  - e2)(l. 0051353. - 0.0025731sm2<p) 


a = 6 378  249.20  m,  ez  = 0.006  803  4877,  e'2  = 


1 — e2 

1.  Montrer  que  les  coordonnées  du  point  A sont  : X = 157  833.48  m , Y = 4 078  512.97  m, 
on  justifie  les  résultats. 

2.  Soit  le  point  B de  coordonnées  (X  = 160  595.98m;  Y = 4 078  564.53m).  Sachant 
que  B est  situé  sur  le  même  parallèle  que  A,  calculer  la  longitude  A'  de  B. 

3.  Calculer  le  gisement  G et  la  distance  AB. 

4.  Sachant  que  la  convergence  des  méridiens  7 est  donnée  par  fgq  = (À  — Ào )sinp  et 

qu’on  néglige  le  Dv,  calculer  l’azimut  de  la  direction  AB. 

5.  Calculer  l’azimut  de  B vers  A en  négligeant  le  Dv  de  B vers  A. 

6.  En  calculant  les  coordonnées  UTM  de  A et  R,  on  trouve  respectivement  Xa  = 
657  770.34 m,  YA  = 4076  891.20  m;  XB  = 660  531.74 m,  YB  = 4 076  942.76m,  Calculer 
la  distance  AB  par  les  coordonnées  UTM.  En  déduire  l’erreur  relative  sur  la  distance  en 
utilisant  les  coordonnées  de  l’UTM  tronquée. 
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3.10  Les  Transformations  de  passage  entre  les  Sys- 
tèmes Géodésiques 

Problème  n°  1 : Soient  les  trois  tableaux  ci-dessous  des  coordonnées  3D  respective- 
ment dans  les  systèmes  S'1  et  S 2 et  à transformer  dans  le  système  S 2 : 


Nom . 

X(m) 

Y (m) 

Z(m) 

1 

4300  244.860 

1062  094.681 

4 574  775.629 

2 

4 277  737.502 

1 115  558.251 

4 582  961.996 

3 

4 276  816.431 

1081197.897 

4 591886.356 

4 

4315183.431 

1135  854.241 

4 542  857.520 

5 

4 285  934.717 

1110917.314 

4 576  361.689 

6 

4217271.349 

1193  915.699 

4 618  635.464 

7 

4 292  630.700 

1079  310.256 

4579117.105 

Nom . 

X(m) 

Y (m) 

Z 

(m) 

1 

4300  245.018 

1062  094.592 

4 574 

775.510 

2 

4 277  737.661 

1 115  558.164 

4 582 

961.878 

3 

4 276  816.590 

1081197.809 

4 591 

886.238 

4 

4315183.590 

1135  854.153 

4 542 

857.402 

5 

4 285  934.876 

1110917.227 

4 576 

361.571 

6 

4217271.512 

1193915.612 

4 618 

635.348 

7 

4 292  630.858 

1079  310.168 

4 579 

116.986 

Nom . 

X(m) 

Y (m) 

Z(m) 

A 

4351694.594 

1056  274.819 

4 526  994.706 

B 

4319  956.455 

1095  408.043 

4 548  544.867 

C 

4303467.472 

1110  727.257 

4 560  823.460 

D 

4 202  413.995 

1221146.648 

4 625  014.614 

1.  Déterminer  les  paramètres  du  modèle  de  Bursa-Wolf  à 7 paramètres. 


2.  Calculer  les  coordonnées  3D  des  points  du  troisième  tableau  dans  le  système  S 2. 
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CHAPITRE  3.  EXERCICES  ET  PROBLEMES 


3.11  Notions  sur  le  Mouvement  d’un  Satellite  Ar- 
tificiel de  la  Terre 

Exercice  n°  1 : 1.  Montrer  que  : r = a(l  — ecosE). 

2.  Démontrer  à partir  des  formules  du  cours  la  relation  : 

v 11  + e E 

tg2  ~ V 1 - 2t9~2 

Aide  : exprimer  tg(v/2)  en  fonction  de  tgv. 

Exercice  n° 2 : A partir  de  l’expression  de  Xc,  monter  que  Xc  vérifie  l’équation  du 
mouvement  non  perturbé  pour  la  composante  X,  soit  : 

Ac  + ^XC  = 0 


Problème  n°l  : La  Terre  est  supposée  sphérique,  homogène  de  rayon  R = 6 371  000  m. 
Le  produit  de  la  constante  universelle  de  gravitation  terrestre  G par  la  masse  M de  la 
Terre  soit  GAI  = 3.986  005  1014  m3,V2.  Un  satellite  géodésique  a une  trajectoire  telle 
que  son  altitude  maximale  est  1100  km  et  son  altitude  minimale  800  km. 

1.  Donner  la  période  de  ce  satellite. 

2.  Quelle  est  l’excentricité  de  sa  trajectoire? 

3.  On  mesure  la  distance  du  satellite  à une  station  au  sol  de  latitude  43°,  5 et  d’altitude 
nulle,  lors  du  passage  du  satellite  à la  verticale  de  la  station,  soit  D = 812  000  m. 

a - Quelle  est  l’anomalie  vraie  du  satellite  à cet  instant,  sachant  qu’il  vient  de  passer 
au  périgée. 

b - Combien  de  temps  s’est  écoulé  depuis  le  passage  au  périgée  ? 

Problème  n° 2 : Une  comète  décrit  autour  du  Soleil  une  ellipse  d’excentricité  e de 
demi-grand  axe  a et  de  demi-petit  axe  b où  le  Soleil  occupe  un  des  foyers.  L’équation  de 
l’orbite  de  la  comète  en  coordonnées  polaires  est  donnée  par  : 

a(l  — e2) 

r = 

1 + ecosv 

avec  r la  distance  Soleil-  comète. 

1.  Déterminer  les  distances  ta  et  rp  lorsque  la  comète  est  à l’apogée  et  au  périgée  en 


fonction  de  a et  e. 
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2.  La  comète  de  Halley  a une  orbite  fortement  excentrique  : son  apogée  est  à 0.53 
U A du  Soleil  et  sa  périgée  est  à 35.1  U A.  Calculer  e. 

3.  En  utilisant  la  loi  des  aires  et  la  troisième  loi  de  Kepler,  montrer  que  la  constante 
des  aires  C est  exprimée  par  : 

p 62 

C2  = —GM 
a 

où  G,  M désignent  respectivement  la  constante  de  la  gravitation  universelle  et  la  masse 
du  Soleil. 

4.  On  pose  : u = Donner  l’expression  du  carré  de  la  vitesse  v2  de  la  comète  en 

r 

fonction  de  u et  — . Montrer  que  v 2 peut  s’écrire  sous  la  forme  : 
dv 

v2  = GM  (--- 
\r  a 

5.  Déterminer  l’expression  du  rapport  des  vitesses  à l’apogée  et  au  périgée  — en 

vP 

fonction  de  e. 

6.  Calculer  numériquement  ce  rapport  pour  le  cas  de  la  comète  de  Halley. 

On  donne  : 

- 1 U A = 149  597  870  km, 

- G = 6.672  x KT11  m3.kg~1.s~2, 

- M = 1.9891  x 1030%. 


CHAPITRE  3.  EXERCICES  ET  PROBLEMES 


Troisième  partie 


Théorie  des  Erreurs 
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Chapitre 


Exercices  et  Problèmes 

4.1  Exercices  et  Problèmes 

Exercice  n°  1 : Soient  U un  ouvert  convexe  d’un  espace  de  Banach[|]  réel  E c’est- 
à-dire  un  espace  vectoriel  normé  complet  sur  K et  / une  fonction  à valeurs  réelles, 
différentiable  et  convexe  dans  U.  Montrer  que  si  f'{x o)  = 0 en  un  point  xq  G U , alors  / 
a un  minimum  absolu  en  xq. 

Exercice  n° 2 : Montrer  que  dans  un  espace  de  Banach  réel  E,  la  fonction  / = ||  .||2 
est  strictement  convexe,  c’est-à-dire,  Va  £]0, 1[,  f(ax  + (1  — ot)y)  < af(x)  + (l  — a)f(y), 
pour  tout  couple  (x,y)  G E2. 

Aide  : utiliser  l’identité  remarquable  : 

\\olx  + (1  - a)y||2  = a||x||2  + (1  - a)||y||2  - a(l  - a)\\x  - y |[2 


Exercice  n° 3 : On  note  F une  surface  de  M3  définie  par  la  représentation  paramé- 
trique : 

OM  = (ai(u,v),a2(u,v),a5(u,v))T 

où  u et  v sont  deux  paramètres  réels.  On  se  donne  un  point  P(x,  y , z)  G M3. 

1.  Donner  une  condition  géométrique  portant  sur  le  plan  tangent  à E au  point 
Mo(uo,  uo)  pour  que  la  différentielle  de  la  fonction  (it,  v)  — > <p(u,  v)  = ||  OP  — OM (u,  v)  || 2 
soit  nulle  en  Mq(uq,vo). 


1.  Stefan  Banach  (1892-1945)  : mathématicien  polonais. 
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CHAPITRE  4.  EXERCICES  ET  PROBLEMES 


Exercice  n° 4 : Soient  U un  ouvert  convexe  d’un  espace  de  Banach  réel  E et  f une 
application  différentielle  de  U dans  M. 

1.  Montrer  que  / est  convexe  dans  U si  et  seulement  si  : 

f(x)  > f(x0)  + f'(x 0){x  - Xq) 


pour  tout  couple  de  points  x,  xo  G U. 

2.  On  suppose  E = Mn  et  / de  classe  C 2 soit  deux  fois  différentiable  et  /”  continue  ; 
pour  x G U,  soit  ipx  la  forme  quadratique  définie  par  : 


Tx{h)  = ^2 

i,j= 1 


d2f 

dxidxj 


( x)hihj , 


h = (hi,h2,...,hn)  G Mn 


Montrer  que  / est  convexe  dans  U si  et  seulement  si  (px  est  positive  pour  tout  x G U 
soit  (px(h)  > 0 pour  x € U et  h € Rn. 

Exercice  n° 5 : Soit  un  triangle  ABC , on  observe  les  angles  A,  B,  C et  les  côtés 
BC  = a,  AC  = b et  AB  = c : 

Â = 43.7716  0gr  a ^ = 3.1  drngr 
Ê = 98.3904  3 gr  = 3.1  dmgr 

C = 57.8385  8 gr  = 3.1  dmgr 

a = 333.841  m,  aa  = 0.005  ru. 
b = 525.847  m,  <rb  = 0.010  m 
c = 414.815  m,  <rc  = 0.005  m 


1.  Calculer  les  angles  et  les  côtés  compensés. 

2.  Calculer  les  poids  de  l’angle  A et  du  côté  a. 

3.  Déterminer  une  estimation  du  facteur  de  variance  unitaire. 


Problème  n°  1 : Les  directions  suivantes  sont  observées  respectivement  aux  stations 
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A,  B,  C et  D d’un  quadrilatère  ABDC  comme  suit  : 


Station  A = 


vers  B : 0.0000  0 gr 
vers  C : 74.1666  7 g?’ 


Station  B 


vers  D : 0.0000  0 gr 
< vers  C : 82.4608  0 gr 
vers  A : 170.6253  lgr 


Station  C 


vers  A:  0.0000  Ogr 
< vers  B : 37.6709  9 gr 
vers  D : 85.0830  2 gr 


Station  D = 


vers  C : 0.0000  Ogr 
vers  B : 70.1280  9gr 


Les  observations  sont  non  corrélées,  l’écart  quadratique  moyen  de  ces  observations  est 
identique  et  vaut  ad  = 6.2  dmgr. 

1.  Compenser  les  directions  et  calculer  leurs  poids  et  celui  de  l’angle  CB  A. 

s 2 

2.  Calculer  l’estimateur  s 2 du  facteur  de  variance  unitaire  et  celui  de  — r. 

az 

3.  Des  observations  de  nivellement  ont  été  effectuées  sur  les  lignes  ABC  et  BCD. 
Les  différences  d’altitudes  observées  sont  les  suivantes  : 


Ha  — Hb  = 0.509??r 
Hb  — Hb  = 1.058  m 
Ha  — Hq  = 3.362  m 
Hb  — Hq  = 1.783  m 
HB-  Hc  = 2.829m 

Les  observations  sont  non  corrélées  et  de  précision  identique.  Compenser  les  observations 
ci-dessus  et  calculer  un  estimateur  du  facteur  de  variance  unitaire. 


Problème  n° 2 : 1.  Montrer  que  dans  un  cheminement  altimétrique  de  précision,  le 
poids  de  l’observation  entre  deux  repères  est  inversement  proportionnel  de  leur  distance 
en  supposant  l’égalité  des  portées  et  que  les  observations  sont  non  corrélées. 


2.  Une  polygonale  ABCD  (voir  Fig.  4.1.1  ) a été  observée  par  le  nivellement  de  pré- 


cision. L’instrument  utilisé  a une  précision  de  2 mm  par  km.  Les  observations  considérées 
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C 


B 


Figure  4.1.1  - La  polygonale  observée 
non  corrélées  sont  les  suivantes  : 

Hq  — Ha  = 1.878  m,  AC  = 6.44  km 
Ho  — Ha  = 3.831  m,  AD  = 3.22  km 
Ho  — Hq  = 1.954  m,  CD  = 3.22  km 
H b — Ha  = 0.332  m,  AB  = 6.44  km 
Ho  — Hb  = 3.530  m,  BD  = 3.22  km 
Hc  — H b = 1.545  m,  BC  = 6.44  km 

L’altitude  du  repère  A est  de  3.048  m et  non  entachée  d’erreurs.  Calculer  par  com- 
pensation des  observations  les  altitudes  des  repères  B,  C et  D et  leurs  écarts-types. 

3.  Calculer  l’écart-type  de  la  différence  d’altitudes  entre  les  repères  C et  D. 

4.  Donner  une  estimation  de  la  précision  par  km  du  nivellement  effectué. 

Problème  n°3  : On  veut  étalonner  un  anéroide,  appareil  donnant  la  pression  de 
l’air,  par  la  formule  : 

D = d + at  + 7 

où  a,  7 sont  deux  constantes,  t la  température  en  degrés  centigrades.  Les  paramètres  d 
et  D sont  lus  respectivement  de  l’anéroide  et  à partir  d’un  baromètre  en  mercure,  et 
exprimés  en  mm  Hg. 

Pour  déterminer  a et  7,  des  lectures  ont  été  prises  à différentes  températures  (voir 
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tableau  4.1.1  ).  Ces  observations  sont  non  corrélées.  L’écart-type  de  la  lecture  de  d est  de 


t 

d 

D 

0 Centigrade 

mm.  H g 

mm  H g 

6.0 

761.3 

762.3 

10.0 

759.1 

759.5 

14.0 

758.4 

758.7 

18.0 

763.1 

763.0 

Table  4.1.1  - Table  des  observations 


0.14  mm  H g]  t et  D sont  supposées  sans  erreurs. 

1.  Calculer  par  la  méthode  des  moindres  carrés  les  constantes  a et  7. 

2.  Estimer  le  facteur  de  variance  unitaire. 

3.  Déterminer  la  variance  et  la  covariance  de  a et  7. 


Problème  n°4  : En  statistiques,  la  loi  normale  est  une  famille  de  distributions  de 
probabilités  caractérisées  par  la  fonction  de  densité  : 

1 

p(x,n,a)  = e 2^ 

v27Tcr 

où  fi  est  la  moyenne  et  cr1 2  la  variance.  On  note  par  l(x,fi,cr)  = Logp(x,  fi,cr),  soit  : 


l(x,  fi,  a)  = —Loger  — 


(x  - fif 
2cr2 


Soit  X une  variable  aléatoire  ayant  comme  fonction  de  densité  p{x,  fi,  cr).  On  rappelle  les 
opérateurs  suivants  espérance  mathématique  ou  moyenne  et  variance  : 

r*+oo 


/+00 

f(x)p(x,  fi,  o)dx 

-OO 

V(f(X))  = E[(E[f(X)}  - f(X))2] 


On  donne  la  formule 


r+ 00 

:/o 


e~u  du  = 


7 r 


1.  Montrer  que  : 

/+00 

p(x,  fi,  o)dx  = g 

-00 

/ + OO 

( x — fi)2p(x,  fi,  a)da 

-00 


Ix  = a 
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2.  Montrer  que  : 


3.  Calculer 


dl_  dl_ 
dg  Ôa' 


Sy/ÎT 


4 


dl  dl 

4.  On  pose  6 = (a,  a).  Soit  Tg  l’espace  engendré  par  , — — ).  On  définit  sur  Tg  l’opé- 

dfi  da 


rateur  Tg  X Tg 


à A,  B deux  variables  aléatoires  G Tg  : 


<A,B  >=  E[A(x)B(x)] 


Justifier  qu’on  peut  écrire  : 

E[A(x)B(x)]  = Cov(A(x),  B(x))  = E[{E[A(x)}  - A(x))(E[B(x)\  - B(x))] 
5.  Montrer  que  < > définit  un  produit  scalaire  sur  Tg. 


6.  On  pose  : e\ 


— et  eo  = On  définit  le  tenseur  métrique  sur  Tg  par  : 
dfj.  da 

9ij  — ^ Ci)  Gj  ^ 


Montrer  que  la  matrice  g = (gij)  est  donnée  par  : 


et  que  la  première  forme  fondamentale  sur  Tg  s’écrit  : 


ds2 


u 


(dfi2  + 2 da2) 


Problème  n° 5 : Soit  un  triangle  de  côtés  a,  b,  c et  d’angles  A.  B et  C.  On  se  propose  : 
- d’estimer  à,  b et  c,  et  les  variances  de  ces  déterminations.  Les  observations  sont  : 

a = 96.48  mm 
b = 115.50  mm 

< A = 63.042  gr  (4.1.1) 

B = 99.802  gr 
C = 37.008  gr 


On  choisit  ici  comme  unités  normalisées  le  décimillimètre  (0.1  mm)  pour  les  mesures  de 
distances,  et  le  décimilligrade  (0.1  gr)  pour  les  angles. 
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On  prend  les  poids  égaux  aux  inverses  des  carrés  des  emq  de  chaque  observation.  On 
donne  la  matrice  des  poids  P : 


P = 


( 0.277 
0 
0 
0 

V 0 


0 

0.160 

0 

0 

0 


0 0 

0 0 

1.524  0 

0 1.524 

0 0 


0 ^ 
0 
0 
0 

1.524/ 


On  prendra  comme  valeurs  approchées  des  inconnues  a o = a ; 6q 


b; 


cq  = a 


sinC 

sinA 


1.  Ecrire  les  paramètres  observées  et  les  valeurs  observées  des  inconnues  dans  les 
nouvelles  unités. 

2.  Soit  X = ( a,b,c ) le  vecteur  des  inconnues.  On  adopte  le  système  suivant  liant  les 
inconnues  aux  observables  : 


a = a 


6 = 6 


b2  + (?  — à2  ■ 

Arccos ; = A 

26c 


c2  + a2  - b2  ■ 

Arccos = n 

2 cà 


(4.1.2) 


Arccos 


à2  + b2  - c2 


2àb 


= C 


Ceci  étant,  on  posera  pour  les  grandeurs  à déterminer  : 


à = ao  + da  = a + da 
b = 6o  + db  = b + db 
è = cq  + de 


à = a + va 
b = b + vb 
À = A + va 
B = B + vb 
C = C + vc 


et  pour  les  grandeurs  observées  : 
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En  linéarisant  la  troisième  équation  de  (4.1.2),  montrer  que  l’équation  d’observation 
s’écrit  : 


1 ao  2000  1 ai  + bi  — ci  2000  1 ai  + ci  — bi  2000 

da : — 7— — -ré db— 


sinA  boco  it  sinA  26qCo  tt  sinA  26oCq  n 
où  : 

^ _ 6q  + Cg  - ag  - 2b0cocosA 
2bç,cosinA 

(étant  entendu  qu’on  exprime  va  en  dcgr). 


, , 2000 
de  = —ka \~va 


TT 


3.  Montrer  que  le  système  des  moindres  carrés  AX  = L + V s’écrit  : 


( L 

0. 

0.  N 

( °'  1 

^ Va) 

0. 

1. 

0. 

^ dti' 

0. 

Vb 

1.00375 

-0.83924 

0.00143 

db 

= 

0.97981 

+ 

va 

-1.00571 

1.20285 

-0.66128 

\dcj 

-2.88449 

VB 

\ 0.00094 

-0.36239 

0.65918  ) 

\ 0.42396  ) 

\vc) 

4.  Résoudre  le  système  précédent  par  la  méthode  des  moindres  carrés  et  montrer  que  la 
matrice  normale  N = À1  PA  est  donnée  par  : 


N = 


( 3.35605 


V - 


-3.13044  1.01750  \ 

3.64132  -1.57937 

1.32971 


5.  Montrer  que  : 


X = 


/ +0.6297l\ 
-0.90962 
\ 0.94782  J 


6.  Déterminer  les  variances  des  inconnues  cr^,  (i‘l  et  aj?. 


Exercice  n°l  : On  considère  (u,v)  G M2  et  on  définit  la  fonction  par  : 

f[u , v)  = tt4  + 6 uv  + 1.5v2  + 36r>  + 405 

1.  Chercher  les  points  critiques  réels  de  /. 

2.  Montrer  que  le  point  x*  = (u,v)  = (3,  —18)  est  un  point  minimum  de  /. 

3.  Montrer  que  le  Hessien  de  / est  une  matrice  définie  positive  si  u2  > 1 et  indéfinie 
si  u2  < 1. 
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4.  Montrer  que  la  formule  de  récurrence  de  Newton  s’écrit  avec  J = 1.5(it?  — 1)  : 


ul 


'U’k+l  — 


J 


Vk+l  — 


2u\  + 18  u\ 

J 


Problème  n°6  : Soient  le  plan  (P)  et  la  sphère  (S2)  d’équations  respectivement  : 
x + y + z = 1 et  x2  + y2  + z2  = 1.  On  veut  chercher  le  point  M G (S2)  tel  que  sa  distance 
au  plan  (P)  soit  maximale. 

1.  Montrer  que  la  distance  d’un  point  M (X,Y,  Z)  au  plan  ( P ) est  donnée  par  : 

d=  \X  + Y + Z 


2.  Pour  répondre  à la  question  posée  ci-dessus,  on  considère  la  fonction  : E(x,  y,  z,  A)  = 
— (x  + y + z — l)2  — À(.x2  + y2  + z2  — 1).  Ecrire  le  système  d’équations  donnant  les  points 
critiques  de  E qu’on  note  par  (1). 

3.  Montrer  que  si  A = — 1,  on  arrive  à une  contradiction.  On  suppose  que  À 7^  — 1. 
Que  représente  le  cas  À = 0. 

4.  On  suppose  que  A ^ {—3,  —1,  0}.  Résoudre  le  système  (1).  Soit  le  point  M2  tel  que 
ses  coordonnées  sont  négatives. 

5.  Montrer  que  la  matrice  hessienne  de  E pour  s’écrit  sous  la  forme  : 


/V 

-2 

-2> 

H = 

-2 

-2 

V-2 

-2 

M2/ 

avec  n=  1 + V3 


6.  Si  on  pose  U = (X,  Y,  Z)T  G (S2).  Montrer  que  UT .H.U  = 2 [3  + ^3  -(X  + Y + Z)2} . 
En  déduire  que  UT .H.U  > 0 pour  tout  U / 0 G (S2). 

7.  Montrer  que  pour  le  point  M2,  on  obtient  un  minimum  strict  de  E.  A-t-on  répondu 
à la  question  du  problème. 


Problème  n°7  : Dans  le  plan  affine  V , on  a mesuré  trois  distances  planes  entre 
un  point  inconnu  P(Xi,  X2)  vers  trois  points  connus  13  dans  trois  directions 

différentes.  On  considère  le  modèle  non  linéaire  de  Gauss-Markov  défini  par  : 

c (X)  = L-e,  e G Af(0,  T) 


avec  : 


- L : le  vecteur  des  observations  (3xl)  = (L\,  L21L^)T  ; 
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- X : le  vecteur  des  inconnues  (2xl)  = (X\,  X2)T  ; 

- e : le  vecteur  des  erreurs  (3x1)  = (ei,e2,e3)r  suit  la  loi  normale  A/"(0, T)  avec 
E(e)  = 0 et  T = E(eeT ) la  matrice  de  dispersion  ou  variance,  on  prendra  T = (Tq.P”1, 
P est  la  matrice  des  poids  égale  à la  matrice  unité  1%,  ao  une  constante  positive  ; 

- ()  : est  une  fonction  donnée  injective  d’un  ouvert  U C R2  — > R3  définie  par  : 


C(X)  = C(*1,*2)  = 


A [(Xi-ai)2  + pf2-6i)2]\ 

2 [(*1  - û2)2  + (X2  - b2)2] 
\l[(Xl-a3)2  + (X2-btf]) 


On  prendra  comme  composante  Lj  du  vecteur  observation  la  quantité  Lj  = 


D 2 


1 observée 


dÇ  d( 

1.  Montrer  que  les  vecteurs  , sont  linéairement  indépendants  en  chaque 


point  X G U. 

2.  Montrer  que  les  fonctions 

3.  Posons  : J = \\L  — £(X)||2 


dX1  ’ dx2 
d2( 


dXidXj 


sont  continues  sur  U pour  i , j G {1,2}. 


Calculer  les  coefficients  de  la  matrice  ( 
4.  Soit  la  matrice  carrée  définie  par  : 


d2J 


dXidXj 


1,2}. 


, d((X)  dÇ(X) 

ff(x)  = ( ffij ) avec  gi:j  =<  — , > 


dX,  ’ dx. 


Calculer  les  coefficients  gtj. 

5.  Introduisons  la  matrice  B définie  par  : 


i = 1,2 
J = 1,2 


B(X,  L)  = ( Bij ) avec  BV]  = gij-  < L - ((X), 


d2( 


dXidXj 


> 


i = 1,2 
J = 1,2 


Calculer  les  éléments  de  la  matrice  B et  montrer  qu’elle  est  définie  positive. 
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